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We present some analytical results for the stochastic sandpile model studied earlier by 
Manna. In this model, the operators corresponding to particle addition at different sites 
commute. The eigenvalues of operators satisfy a system of coupled polynomial equations. 
For an LxL square, we construct a nontrivial toppling invariant, and hence a ladder operator 
which acting on eigenvectors of the evolution operator gives new eigenvectors with different 

'^ . eigenvalues. For periodic boundary conditions in one direction, one more toppling invariant 

can be constructed. We show that there are many forbidden subconfigurations, and only an 
O [ exponentially small fraction of all stable configurations are recurrent. We obtain rigorous 

lower and upper bounds for the minimum number of particles in a recurrent configuration, 

^ ■ and conjecture a formula for its exact value for finite-size rectangles. 

[£ : PACS nos: 05.20.Ln, 02.50.Ga, 02.10.Pk, 81.05.Rm 
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§ : 1 Introduction 

^ ■ Many cellular automaton 'sandpile' models have been studied in recent years as simple the- 

c^ ■ oretical models of self-organized criticality. The prototype of these models is the Abelian 

sandpile model (ASM), proposed by Bak et al in 1987 [jl], 0]. Many properties of this model 
__j , can be determined analytically using the fact that the operators corresponding particle ad- 

C3 I dition in this model generate an abelian group, and the model is equivalent to the q ^ 

^ ■ limit of the Potts model 0. Several variations of the ASM have been studied in the past 

with a view to understand the parameters that determine the different universality classes 



X 



of self-organized critical behavior |^, |^, Q. These include models in which particle transfer 
is directed 0, or models in which the toppling condition or the number of sandgrains trans- 
c^ I ferred depends on the local slope rather than local height. In this respect, it has been realized 

that stochasticity in toppling rules can lead to different critical behavior than models with 
deterministic toppling rules ^ ^ |T0| . 

As a physical motivation for the study of such models, we note that stochasticity in the 
toppling rules may be taken as a simple phenomenological attempt to take into account in 
a theoretical description the variation in shape and smoothness of different grains in real 



granular media. It has been argued that a stochastic one-dimensional sandpile model seems 



to describe well the behavior of experimental results on avalanches in ricepiles |TT| . 

In this paper, we study a class of stochastic sandpile models, the prototype of which 
is the so-called Manna model [§. It was recently proved that this model also shows the 
abelian property, and is a special case of the more general Abelian Distibuted Processors 
(ADP) model 0. In the following, we shall use the terms deterministic ASM (DASM) and 
stochastic ASM (SASM) if we need to distinguish between these classes of models. We will 
show that some of the analytical techniques of DASM are also useful in studying the SASM. 

The plan of this paper is as follows: In section 2, we recapitulate the definition of the 
Manna model. We study this model in two dimensions on a square lattice, but most of 
the treatment can be easily generalized to higher dimensions. In section 3, we describe 
the abelian ring structure of the algebra generated by the particle addition operators. The 
simultaneous eigenvalues of these operators satisfy a system of coupled quadratic equations. 
In section 4, we show that for the case of an L x M rectangle, if L = M, there is a 
nontrivial toppling invariant. We use this to construct a ladder operator that gives (L + 1) 
simultaneous eigenvectors of the particle addition operators from one such eigenvector. For 
periodic boundary conditions in the y-direction, and M = L + 1, we can construct one more 
toppling invariant, and this raises the number of eigenvectors generated from one vector to 
M^. In section 5, we prove the existence of forbidden sub-configurations (FSC's) in the 
model. We describe an algorithm to construct the set of minimal FSC's. In section 6, we 
study the recurrent configurations of the model which have minimum number of sandgrains. 
We prove upper and lower bounds on this number for a finite rectangle, and conjecture that 
the upper bound always coincides with the exact value. 

2 Definition of the Model 

We consider a two-dimensional square lattice. The sites of the lattice will be labelled by 
their integer Euclidean coordinates {x,y), 1 < x < L, 1 < y < M. At each site {x,y), 
there is a non-negative integer Zr^^y called the height of the pile. In addition, at each site, we 
assume that we have a pseudo-random- number generator (PRNG), which gives a number 
lying between and 1 on each request using a deterministic algorithm. The n-th output 
from a PRNG depends only on its initial setting (seed), and n, and not in any way on 
the evolution of the system, or status of other PRNG's. Equivalently, we can think of the 
PRNG's as an infinite stack of independent random numbers at each site. 

If the height at any site does not exceed 1, the configuration is said to be stable. If the 
system is in a stable configuration, we choose a site at random and add a grain of sand 
there, as a result of the addition, height at that site increases by 1. If this makes the site 
unstable, we relax the system using the relaxation rule given below. Once all unstable sites 
are relaxed, we add another particle. We shall assume that the probability of addition is 
nonzero for all sites. It follows from the general theory of ASM's that the steady state of 
the system is then independent of the precise choice of these probabilities. 

The relaxation rule is as follows: Any site whose height exceeds 1, relaxes by toppling, 
and throwing two grains out. Where these grains go is decided by drawing a new random 
number from the local PRNG. In terms of stacks, we pop up a number from the local stack. 



If the number is less that 1/2, one particle is transferred to each of the two horizontal 
neighbors of the site. Else, one particle is transferred to each of the two vertical neighbors. 
If the transfer renders some other sites unstable, these are toppled in turn, until a stable 
configuration is reached. Sandgrains can leave the system if a toppling occurs at a boundary 
site. 

3 The Abelian Ring of Operators 

This model as a special case of the more general Abelian Distributed Processors Model 
(ADP) 1^, and the abelian character of the model becomes obvious. We define particle 
addition operators A^^y which acting on a stable configurations C, gives the stable configu- 
ration obtained after adding a particle at {x,y), and relaxing the system. Here the stable 
configurations are specified by the heights {z^^y}, and the state of the PRNG at each site 
(equivalently, a stack of random numbers at each site). It is easy to see that if we start with 
a configuration with two or more unstable sites, we get the same configuration, independent 
of the order in which the unstable sites are relaxed. Using this property repeatedly, it follows 
that 

[A^,y, Aa;/y] = 0, for all sites {x,y) and {x',y'). (1) 

We may now assume that the internal state of these PRNG's as inaccessible to the 
outside observer, and that there are no observable correlations between different outputs of 
the PRNG's ( i.e. they are 'good' PRNG's). Then the deterministic evolution of avalanches 
in the (sandpile heights+ PRNGs) system may be described equally well as a Markovian 
relaxation of the subsystem described in terms of sandpile heights alone. 

Given a particular state of the PRNG's at each site, addition of two particles at different 
sites gives same final stable configuration, independent of the order of topplings. A different 
state of PRNG's can lead a different final state, but still independent of order of topplings. 
If we do not know the state of the PRNG's, we have to average over their different possible 
states. Then, the final state is not fixed, only the probabilities of different outcomes can be 
given. 

In Manna's original formulation of the model, the transfer direction was decided randomly 
at each toppling. The advantage of the specific implementation used here with a different 
PRNG at each site makes the fact that abelian character of the model obvious. The abelian 
property, once established, is clearly independent of the details of implementation of the 
algorithm. So, one could as well use the same PRNG to decide the direction of particle 
transfer in topplings at different sites. 

We consider the 2-^*^ -dimensional vector space V, whose basis vectors \Z > are labelled 
by the stable height configurations Z = {z^^y}- If the probability that at time t the config- 
uration of the sandpile is Z is denoted by Prob{Z,t), then the state of the system can be 
represented by a vector \P{t) >G V 

\P{t)>=Y,Prob{Z,t)\Z> (2) 

z 



We define particle addition operators a^ ,, as linear operators acting on V as follows: for all 
configurations 2, a^. j^|2 > gives the probability vector of the system obtained by taking the 
system in the stable configuration Z^ adding a particle at site (x, y\ and relaxing the system 
until a stable configuration is reached. For stochastic toppling rules, the resulting state is not 
necessarily a basis vector corresponding to a unique configuration, but a linear combination 
of different basis vectors. Note that action of any of these operators on a given configuration 
gives a unique probability state vector, and action of these {a^ y}'s is deterministic, in the 
same way as evolution of the wave-function in quantum mechanics is deterministic. The 
Abelian property of the Manna model means that the probabilities of different final stable 
configurations are independent of the order of topplings and particle additions. 

The commutativity of the operators {A(x, y)} implies that the operators {Si-x,y} also 
commute with each other: 

[dix^y.Six'^yi] = 0, for all sites {x,y) and {x',y'). (3) 

This abelian property simplifies the analysis of the model considerably. There are addi- 
tional relations satisfied by the operators {a^^y}. Adding two particles at a site will cause it 
to topple, whatever the initial configuration, and then with equal probability two particles 
are transferred either horizontally, or vertically. This implies that the operators {a^, j^} satisfy 
the equations 

^l,y = ^{^x-i,yiix+i,y + a.x,y-iax,y+i) , for all X, y, {I < X < L, 1 < y < M) . (4) 
In writing these equations, we have assumed the boundary conditions 

ao,y = SLL+i,y = 1, ioT 1 < y < M (5) 

ax,o = a^,M+i = 1, for 1 < a: < L (6) 

If we assume periodic boundary conditions in the y-direction, the Eqs.(6) are replaced 
by the equations 

Six,y = 8Lx,y+M, for all x, 1 < a; < L, all y. (7) 

Any vector v G V can be expressed as P(v)|<I> >, where -P(v) is a polynomial in {a^ j,}, 
with the highest power of any a^. ^ being 1, and 1$ > is the vector corresponding to config- 
uration with all sites empty. We can multiply two such polynomials Pi and P2, and reduce 
the highest power of operators Six,y using the reduction rules Eq.(4). 

We may think of the operators {six,y} as abstract symbols, and study the algebra of 
polynomials in these variables defined by the usual addition, and multiplication. Using the 
reduction rules (4), the maximum degree of any symbol in the polynomials can be reduced 
to 1. Clearly, the set of all such polynomials forms a commutative ring. 

It is instructive to compare this situation with that of the DASM. In the latter case 
also, there are reduction rules which reduce large powers of abelian operators {a^;,^} to lower 
ones, but the expression on the right-hand side of Eq. (4) consists of a single term. Then 
the operators define a finite-dimensional abelian semi-group, and if we restrict to the set of 



recurrent states, we get an abelian group with an well-defined multiplicative inverse 0. In 
the present case, the reduction rules give us not a single term, but a linear combination of 
terms. To form an algebra closed under addition and multiplication, we have to include all 
linear combinations of terms. Thus, the number of elements of the algebra is infinite. 

The operators {aa;^^} have an obvious representation as 2^^^ x 2^^ matrices ( defined by 
their action on basis vectors of V). These matrices are in general not diagonalizable, only 
reducible to the Jordan canonical form, and there may not be eigenvectors corresponding 
to each repeated eigenvalue. For the DASM, all transient states correspond to the zero 
eigenvalue of the operators {aj. y}, and all other eigenstates are of modulus 1. For the SASM, 
there may be a set of configurations such that if the system is any one of them, with a finite 
probability, it makes a transition within the same set, and with a finite probability it gets out. 
If the system can never get back to these states, these states are transient, but there need 
not be any zero eigenvalue of the operators {aj. y} corresponding to such transient states. 
Also, amongst the recurrent states, one can have two different initial probability vectors that 
yield the same resultant vector. Thus, there can be eigenvectors of zero eigenvalue in the 
recurrent sector, and in general, the inverse operator a~^ need not exist, even if we restrict 
ourselves to the set of recurrent configurations of the model. 

Consider one of the operators, say ai i. Let Oi^i be one of the eigenvalues of ai i. Let 
Vi be the subspace of V spanned by the (right) eigenvectors of ai^i corresponding to the 
eigenvalue oi^i. There is at least one such eigenvector, so Vi is non-null. We pick one of the 
other addition operators, say ai 2. From the fact that ai 2 commutes with ai 1, it immediately 
follows that ai^2 acting on any vector in the subspace Vi leaves it within the same subspace. 
Diagonalizing ai_2 within this subspace, we construct a possibly smaller, but still non-null, 
subspace V2 which is spanned by simultaneous eigenvectors of ai 1 and ai 2 with eigenvalues 
ai^i and ai,2- Repeating this argument with other operators, we can construct vectors which 
are simultaneous eigenvectors of all the {aL^^y}. Let |\I' > be such a vector, with 

a^.,y|^ >= a^^y\^ >, for all {x,y). (8) 

Then the eigenvalues {ax,y} satisfy the coupled set of quadratic equations 

<^l,y = ■7:{0'x-i,yax+i,y + ^x^y-iax^y+i) , for all x, y; l<x<L,l<y<M. (9) 

These are LM quadratic equations in the LM complex variables {a^-^y}. In general, they 
have 2^^ solutions. It seems difficult to write down the general solution of these equations 
explicitly. However, it is easy to see that all ax^y = 1 is certainly a solution of these equations. 
The corresponding eigenvector |\1' > gives the steady state. Unlike the DASM case, here we 
are not able to write down explicitly the probabilities of different recurrent configurations in 
the steady state. 

4 Toppling Invariants 

Toppling invariants have been found to be very useful in characterizing the structure of 
recurrent states in the DASM case [0. Consider a DASM on N sites, whose configuration is 



specified by the heights {zi}, i = 1,2 . . . N, and the topphng matrix is A. For any arbitrary 
configuration Z = {zi}, a topphng invariant I is defined as a integer function hnear in the 
height variables by the equation 

I{Z) = J2giZi, {mod d). (10) 

i 

where d is some integer, and gi are integers chosen so that / does not change under topphng 
at any of the sites j. The N different choices of j give N conditions to be satisfied by 
the as yet undetermined coefficients {gi}- For suitable choice of d, these equations have a 
non-trivial solution, and in fact one can find a minimal set of toppling invariants in terms of 
the Smith normal form of the matrix A [l^ . 

For the present model, we define a toppling invariant /i by the equation 



h{Z)=^g{x,y) z^^y, {mod d). (11) 

x,y 

For this to be unchanged if there is a toppling at (x, y), we must have 

'^g{x,y) = g{x-l,y) + g{x + l,y), {mod d). (12) 

and 

2fl'(a;,2/) = fi'(x,y- 1) +5f(x,|/ + l), {mod d). (13) 

These are 2LM equations for the LM + 1 unknown coefficients g{x, y) and d. In general, we 
would not expect any nontrivial solution. Interestingly, for the special case L = M, these 
equations have a simple nontrivial solution ioi d = L + 1: 

g{x,y) = xy (14) 

We show below that toppling invariants can also be used as ladder operators, which 
acting on an eigenvector of the evolution operator, gives rise to another eigenvector of the 
same with a different eigenvalue. For the DASM, this is not very important, as the evolution 
operator can be diagonalized directly, and all the eigenvectors are easily determined. For 
the SASM, the eigenvectors are not easily determined, and the fact that toppling invariants 
can be used to construct ladder operators provides useful information. 

We define an operator S, diagonal in the configuration basis by the equation 

S = exp{ j^^-^ ) (15) 

Clearly, S is also invariant under toppling. Simple algebra then gives 

1 2mxy 

S a^,y S = exp{——) a^^y (16) 

Let |\E' > is a simultaneous eigenvector of the operators {sL^^y}, with eigenvalues {a^-y}, 
1 < X < L,l < y < M. Then, it follows that S|\E' > is also a simultaneous eigenvector of all 
these operators with eigenvalues {a^.y}, where 

6 



2Tdxy 
ax,y = a^,y exp[j-^\. (17) 

If we use periodic boundary conditions in the y-direction, and L = M — 1, we can define 
another topphng invariant I2 

I2{Z) =Y,xz^^y, modM. (18) 

x,y 

Thus, for a lattice of size L x (L + 1), we have two independent topphng invariants Ji 
and J2. More generaUy, if M 7^ L + 1, and the greatest common divisor of M and L + 1 is 
R, then both Ji and I2 are invariant modulo R. 

Using these symmtries together, we can construct a set of M^ solutions of the coupled 
eigenvalue equations (9) when M = L + 1, for any solution {a^^y}. These solutions are 
parameterized by two integers k and r. It is easily checked that if {ax,y} is a solution of 
Eqs.(9), then so is {ctx^y}, where 

ax,y = ax,yexp[27ii{k + ry)x/M], < k < M, 0<r < M (19) 

We can group configuration having the same values of Ji, I2 into same equivalence class. 
There are M^ such equivalence classes. Clearly, no value of Ji or I2 has any preferred status, 
and so, in the steady state, each of these equivalence classes must have equal weight. Thus, 
the toppling invariants provide a partial characterization of the probabilities of different 
configuration in the steady state. Also, it is easy to see that all configuration obtained by 
adding a particle at {x, y) to a configuration belongs to the equivalence class of Ji = i,l2 = j, 
would belong to the equivalence class with Ii = i + xy, h = j + x. 

5 Minimal Forbidden Subconfigurations 

Determining the probabilities of different stable configurations in the steady state is not 
easy. However, as a first step towards such a characterization, we show below that most 
of the 2^*^ stable configurations do not occur at all in the steady state, and the number of 
recurrent configurations of the sandpile is only an exponentially small fraction of all stable 
configurations. We do this by showing that there are many local constraints on the allowed 
heights {zx^y} in recurrent configurations. 

The simplest of such constraints is that in any recurrent configuration, any four neighbor- 
ing sites forming a 2 x 2 square cannot be all unoccupied. Thus the local structure shown in 
Fig. la cannot be found anywhere in a recurrent configuration. We shall call such structures 
forbidden subconfigurations (FSC's). In fact, we shall show that an FSC cannot even in an 
intermediate unstable configuration, if we start with a recurrent configuration. 



The proof that Fig. la is an FSC is similar to the proof used in the DASM case [0: We 
start with an initial configuration of the sandpile in which there is no such subconfiguration. 
Now we evolve the sandpile by adding particles at sites, and toppling any sites which become 
unstable. We would like to prove that at no stage, in the intervening stable or unstable 
configurations, would we be able to generate the subconfiguration shown in Fig. la, whatever 
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Figure 1: Some minimal forbidden subconfigurations of the Manna Model 



the sequence of particle additions, and whatever choice we make for each toppling (horizontal 
or vertical transfer). 

To prove, we assume the contrary. Let us treat each particle addition or toppling as a 
seperate relaxation step. Let the FSC appears for the first time after the T-th step. Then, 
after step T — 1, there was no FSC at these sites, and the height at one of these sites must 
have decreased in the T-th step. This implies that the T-th step was a toppling at one of the 
four sites in the 2x2 square. But this toppling, whether transfer was horizontal or vertical, 
would send one particle to one of the neighbors in the square. Since the minimum height 
at the neighbor is zero at time T — 1, its height at the T-th step must be at least 1, which 
contradicts the assumption that the height is at all the four sites at the T-th step. Hence, 
no such T can be found, and the assumption is not possible. 

We can prove that other configurations shown in Fig. 1 are also FSC's similarly. For 
example, to show that sub configuration shown in Fig. lb is an FSC, we again note that if 
this subconfiguration is generated for the first time at the T-th step, then the previous step 
must be a toppling at one of the sites in the FSC. That implies that at the (T — l)-th step, 
the height at one of the neighboring sites was negative, or that the central site was of height 
zero . But the latter case would lead to an configuration shown Fig. la, already proved to 
be a FSC. 

If a subconfiguration cannot be generated on a small lattice, it cannot be generated on 
a bigger lattice either, as the effect of the rest of the lattice on a small part can always be 
reproduced by particles added at the boundaries of the of the part at appropriate points in 
space and time. Thus, every transient configuration of a small lattice is an FSC on a bigger 
lattice. 

Clearly, if F is an FSC, so is any bigger set containing F as a subset. We may therefore, 
without any loss, restrict our discussion to FSC's such that any proper subset of the sites 



does not form an FSC Such FSC's will be called minimal FSC's. 

To form a list of all minimal FSC's, one can start by examining the transient configura- 
tions on small lattices. It is easily seen that all configurations are recurrent if M or L = 1, 
and that Fig. la gives the only transient configuration on a 2 x 2 lattice. One can check 
that going to lattices 2 x M with M > 2 does not generate any new minimal FSC. We then 
enumerate all the transient configurations of the 3x3 lattice. Some of these are disallowed 
by the existence of FSC's of smaller sizes. From the remaining transient configurations, we 
can get a reduced set of minimal FSC's by replacing the FSC by a smaller set with fewer 
sites wherever possible. This gives us a new minimal FSC shown in Fig. lb ( and another 
FSC related to this by reflection symmetry). Going to bigger lattices, we can continue to 
augment our list of FSC's, with larger and larger FSC's. At each stage, we can be sure that 
all minimal FSC's having fewer than a specific number of sites have been enumerated. Some 
of the minimal FSC's constructed in this way are shown in Fig. 1 b-g. If we confined our 
attention to FSC's on stable configurations only, and not insist that they be forbidden in 
unstable states also, these configurations could have been made smaller by omitting all sites 
of height 1 in the figure. 

The existence of FSC's implies local inequality contraints in the height variables that 
must be satisfied at all sites in recurrent states. This implies that the fraction of the 2^^^ 
stable configurations that are recurrent must decrease to zero for large lattices ( exponentially 
fast in the number of sites in the lattice). 

The procedure to construct FSC's is rather tedious, and not very illuminating as the 
number of minimal FSC's is infinite. It seems desirable to have a simple algorithm which 
can identify FSC's. For the DASM, such an algorithm is rather simple: it is known as the 
burning test [|l^, or the script test for some asymmetric graphs 



A simple variation of the burning test seems to work well for a large majority of configu- 
rations: starting with the given configuration, one recursively burns any site i for which the 
height Zi > Min{vi, hi), where Vi and hi are the number of vertical and horizontal unburnt 
neighbors of i respectively. If eventually all sites are burnt, the configuration passes this 
burning test, else not. 

Some not- very-systematic checking on small lattices showed that if a configuration passes 
this modified burning test, it is usually recurrent, and if it fails, it is transient. However, 
there are a few instances where this conclusion is incorrect. For example, the FSC shown in 
Fig. Id burns completely, and hence should have been recurrent according to this test. One 
can also find configurations that do not burn according to this test, but are recurrent. An 
example is given in Fig. 2. It is possible that a small change in the burning condition in this 
test would be make it work without any error. Finding such a refinement remains an open 
problem at present. 

We note that the existence of FSC's is a special feature of these toppling rules. If we had 
chosen the particle transfer to be in the north-and-east or south-and-west directions with 
equal probability, then it is easy to see that starting from any configuration, we can get to 
the configuration with all sites empty, and hence all stable configurations are recurrent. 
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Figure 2: A subconfiguration that does not burn in the modied burning test, but is not an 
FSC. 



6 Minimal Sand Configurations 

Let 7^ be a recurrent configuration of sandpile on a rectangle of size £ x m. We denote 
the number of sandgrains in TZ by N{TZ). The minimum value of N[TZ) over all recurrent 
configurations on the lattice will be denoted by Nmin{£ , fn) ■ In this section we study how 
Nmin{£-,T^) depcuds on a, and m. For large £, m, this should vary linearly with the number 
of sites in the lattice, and we define the minimum possible density of occupied sites in a 
recurrent configuration as Pmin- Clearly, 



Pv 



Limp 



Nmin{£,m) 

im 



(20) 



We have already seen that four adjacent sites forming a 2 x 2 square, if all unoccupied, 
form an FSC, and hence each such group of four sites must have at least one occupied site. 
For large i, m, we can divide the space into non-overlapping 2x2 squares, and so we see 
that 

(21) 



Pr. 



1 
> - 

- 4 



One can get better bounds on Pmin using bigger tiling units. In Fig. 3, we have shown 
the tilings using tiles using 7 and 12 sites each. For the 7-site tile, it is easy to see that we 
get an FSC within the tile, unless the number of occupied sites within the tile is at least 2. 
This implies that pmin > I- For the 12— sited tiling, the minimum number of occupied site 
per tile can easily be checked to be 4. This gives us an improved bound p^nin > 4. 

It is straight-forward to extend this procedure to larger tiles. For any finite tile, there 
are only a finite number of configurations, and by exhaustive enumeration, one can find the 
largest number r, such that all configurations with less than r occupied sites necessarily 
contain an FSC, and are transient. The number of different configurations to be checked 
increases as 2^™", and this brute-force approach does not seem to be very practical for large 
i, m. However, experimenting with rectangular tiles of size i xm leads us to conjecture that 



Nmin{^,m) = f{i,m) 



(22) 
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Figure 3: Tiling of the plane using tiles of 7 and 12 sites. 



where 



f{i,m) 



{i — l){m — 1), ii i or m odd, 
{i'm — i — m + 2), otherwise. 



(23) 



We have checked that the conjecture is true for 1 < £, -m < 4, and if £ < 2 for all m. Note 
that in the limit of large i, m, this implies that 



1 



Pr, 



(24) 



As additional support for this conjecture, we show that Nmin{^, ^ cannot exceed /(£, m). 
This is done by explicitly constructing a recurrent configuration having f{i, m) particles. We 
recall the definition that a configuration is recurrent if it occurs with non-zero probability 
in the steady state of the system. The configuration with all sites occupied can clearly 
be reached from any configuration of the system ( by adding particles only at unoccupied 
sites), and is recurrent. Thus a configuration is recurrent if and only if there is a non-zero 
probability that it is reached in a finite number of steps from the fully occupied configuration. 

We now show that for all values of i and m, starting from a fully occupied configuration 
of size i X m, we can construct a sequence of topplings and particle additions that makes 
the sandpile reach a configuration in which the number of particles is f{i,m). The claim 
is easily checked for i,m < 2. For larger values of i,m, we prove this result by induction. 
Assume that The claim is true for some value {io,mo). We now show that it must then also 
be true for {i = Eq + 2, mg). 

We construct the minimal configuration from the fully occupied configuration of the 
(£o + 2, mo) lattice as follows: In the initial configuration, {zx,y} = 1, for all x,y. We first 
produce a configuration on the rectangle < x < iQ,0 < y < rriQ having only /(^0)'"^o) 
sandgrains. This can be done, by assumption. In this process, topplings at the column 
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X = io would add particles to the column x = io + 1. We do not relax these yet, and add a 
particle to each site of the column x = io + 2, and topple it horizontally. Finally, we relax 
the sites in the column x = io + 1 hj toppling them in the vertical direction only. In the 
X = io + 1 column, we essentially have a one dimensional ASM with deterministic topplings. 
Thus, there can be at most one unoccupied site in the column after relaxation. [If all sites 
in the column are occupied, one vacancy is created on adding one more particle at any site 
in the column, and relaxing, again transferring particles only in the vertical direction]. The 
final number of particles in the configuration is clearly f{£o, mo) + (m — 1) = f{£o + 2, ttiq). 
By symmetry between the horizontal and vertical directions, we have Nmin{£-,f^) = 
N^inim^t). Hence, by induction on (i and m, we conclude that 

Nrmn{£,m)<f{£,m), for all £, m. (25) 

where /(£, m) is the integer function defined by eq(20). A proof of the conjectured equality 
eq.(18) remains an open problem. 

I thank Satya N. Majumdar for a critical reading of the manuscript. 
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